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ABSTRACT 

Context. Mkr 841 is a bright Seyfert 1 galaxy known to harbor a strong soft excess and a variable Ka iron line. 
Aims. It has been observed during 3 different periods (January 2001, January 2005 and July 2005) by XMM for a total 
cumulated exposure time of ~108 ks. We present in this paper a broad band spectral analysis of the complete EPIC-pn 
data sets. 

Methods. We use different methods of data analysis including model-independent methods (spectral ratios, RMS, ...) as 
well as model fitting. We were able to test two different models for the soft excess, a relativistically blurred photoionized 
reflection (ihodel) and a relativistically smeared ionized absorption (ABS model). The continuum is modeled by a simple 
cut-off power law and we also add a neutral reflection. 

Results. These observations conflrm the presence of a soft excess and iron line and reveal extreme and puzzling spectral 
and temporal behaviors. The 0.5-3 keV soft X-ray flux decreases by a factor 3 between 2001 and 2005 and the line shape 
appears to be a mixture of broad and narrow components, the former being variable on small (ks) time scale while 
the later is consistent with being constant. The 2-10 keV spectrum also hardens between 2001 and 2005. We succeed 
in describing this complex broad-band 0.5-10 keV spectral variability using either or ABS to fit the soft excess. Both 
models give statistically equivalent results even including simultaneous BeppoSAX. data up to 200 keV. Both models are 
consistent with the presence of remote reflection characterized by a constant narrow component in the data. However 
they differ in the presence of a broad line component present in but not needed in ABS. Consequently the physical 
interpretation of the line profile variability is quite different, resulting from the variability of the broad line component 
in and from the variability of the absorbing medium in ABS. 

This study also reveals the sporadic presence of relativistically redshifted narrow iron lines, one of them being detected 
at 4.8 keV in the EPIC-pn instruments at more than 98.5% confidence level. If interpreted as the blue horn of a 
relativistically distorted neutral iron line, the large redshift implies the presence of a Kerr black hole. 

Key words. Galaxies: Seyfert; Galaxies: individual: Mkr 841; X-ray: galaxies 



1. Introduction 

Seyfert 1 galaxies emit the bulk of their luminosity in 
the UV and X-ray bands. Their spectral energy distri- 
butions are characterized by two main co mpone nts: a 
UV bump peaking in the UV range (|0'Brien et all Il988t 
iKinnev et al.l . [T99lh . and an apparently non-thermal X-ray 
power law component extending up to at least 100 keV. 
A majority of objects show also the presence of a soft 
X-ray excess with respect to the high energy power la w 
extrapolation (e.g. IWalter fc Fh^ll993t iPage et al.ir2004l ). 

The origin of these components is not well understood. 
The basic paradigm supposes the existence of an accreting 
supermassive black hole. The gravitational energy released 



by the accreting gas is generally thought to be dissipated 
partly in the U V as the r mal h e ating in an opticall y thick 
"cold" plasma (jShieldsl . Il978t iMalkan fc Sargentl . Il982l) 
and partly in X-rays in active blob s of optically thin "hot" 
plasma (e.g. iHaardt et al.l 11994^ . These blobs radiate 
mainly through comptonization of soft photons, possibly 
provided by the optically thick medium. Irradiation of 
the dense gas by hard X-rays produces also an X-ray 
reflection spectrum. The study of this reflected component 
(dominated by the iron Ka emission line and the reflection 
bump peaking between 20-40 keV) appears to be of great 
importance since it has the potential to be a key diagnostic 
of the strong gravity environment of black holes (Fabian 
et al. 2000 and references therein). 

On the other hand, while the origin of the soft excess 
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Table 1. Observation epochs, total duration, % of good exposure time and mean count rates in the EPIC-pn instrument. 



is still not clearly understood, it was recently realized 
that its characteristic temperature (when fitted by e.g. a 
simple black body) is remarkably constant over a wide 
range of AGN luminosities and black hole masses (e.g. 
Czernv ct al. 2003: Gicrli hski fc Dondl2004l : ICrummv et al.l 
2006; Ponti ct al. 2006), favoring an origin through atomic 
processes instead of purely continuum emission. For exam- 
ple, recent studies suggest that an appealing explanation 
could be (photo i onized ) reflection from the accretion disc 
([Crummv et al.l . l2006t ). This spectral decomposition has 
been successful in fitting the XMM-Newton X-ray data of 
many sources like MCG-6-30-15, IH 0707-495, NGC 4051 
or MCG-02-14-009jF abian et al), \2004 iPonti et al.l . l2006t 
iLarsson et al.l l2007l : Porquet, 200^. Moreover, in those 
cases in which a broad Fe line is clearly detected (such as 
in MCG-6-30-15) the model is very attractive because the 
soft excess and broad Fe line are fitted self-consistently 
with the same relativistically blurred reflection model. 
But absorption instead of reflection coul d also repro- 
duce the soft excess (iGierlihski fc Pond. l200l. [200l 
ISobolewska fc Pond . 120071 : ISchurch fc Pond . l2006fl and 
modify the spectral shape close to the iron line to mimic 
the presence of a broad component. Nevertheless the reality 
is almost certainly a co mplex combination of absorption 
and reflection effects (e.g. IChevallier et al.ll2006() and their 
relative importance in the observed spectra has been a 
significant topic of discussion in the recent literature. One 
of the main issues of this debate is the determination of the 
underlying continuum below the iron line in order to per- 
mit a precise measurement of the line broadness, the most 
important signature of the presence of black holes in AGNs. 

Mkr 841 is a bright Seyfert 1 galaxy (z=0.0365, 
iFalco et al.l I1999D . one of the rar e Seyfert Is d e tected 
by OSSE at more than 3 a pohnson et all . Il997t 
IZdziarski et al.l. l2000l). It is kii o wn for its large spe ctral 
variabilitv (|George et all . 119931 : iNandra et all . Il995f) . its 
strong soft excess (t his was the first ob ject where a soft 
excess was observed, lArnaud et al.l I1985D and its variabl e 
iron line (at least on a year time scale. iGeorge iEaD[l99l). 
The later was observed in some cases with a relatively 
large equivalent \ yidth (hereaft e r EW ) of about 400 eV 
(jPav et al.l . Il990l : iBianchi et al.l . l200lD significantly above 
the value pre dicted by standard cold reflection model (e.g. 
IGeorge fc Fabian.il991i) . 



R ecent XMM-Newton observations (jPetrucci et al.l . 

l2002t iLonginotti et al.l . I2004D have revealed a puzzling be- 
havior of the iron line in Mkr 841. Indeed the presence of a 
highly variable, but narrow iron line feature was observed in 
XMM-Newton 2001, completely at odds wit h any currently- 
accep ted interpretation of the line origin (jPetrucci et al.l . 
|2002|) and required further investigation. A re-analysis of 



the XMM-Newton data by ILonginotti et al.l (j2004[ ) pro- 
posed that the line may vary in width rather than in flux. 
Their interpretation then invokes local illumination by a 
flare inducing an hotspot in the inner disc, which then be- 
comes progressively broadened as the disc rotates. On the 
other hand the reflection component, although poorly con- 
strained due to the lack of high signal to noise in the si- 
multaneous BeppoSAX data above 10 keV, was relatively 
large (R>1) confirm ing a previous BeppoSAX observation 
(iBianchi et al.l . 12001*) . Astonishingly, the continuum shape 
and flux kept roughly constant between the two pointings 
as well as during the total (100 ks) BeppoSAX observa- 
tion. For a better understanding of the puzzling spectral 
and temporal behavior of this source, it has been observed 
again in 2005 for a total duration of ~ 75 ks. 

We report here on the detailed spectral analysis of the 
whole XMM-Newton data set including one archive obser- 
vation and the 4 open time pointings done in 2001 and 
2005. 



2. The XMM-iVe?i;ton data 



The first XMM-Newton (|Jansen et al.|[200ll and references 
therein) pointing of Mkr 841 (denoted OBS 1) was done 
on the 13th of January 2001 for ~8 ks as part of the 
guaranteed time program. It was immediately followed 
by the first open time observation. Pue to operational 
contingency, the requested 30 ks were split into two parts, 
on the 13th (OBS 2) and the 14th (OBS 3) of January 2001 
with ~11 and 13 ks duration time respectively. The two 
observations were separated by about 15 hours. The source 
was re-observed 4 years later in January 2005. Here again 
the requested 75 ks were split into two parts due to strong 
proton flares during the observation. About 46 ks of good 
quality data were extracted from this pointing (OBS 4) 
while the 30 ks left were performed in July 2005 (OBS 5). 
Table |T] gives a summary of the XMM-Newton pointings 
with the corresponding dates, observation duration and 
count rates. 

In this paper we (usually) restrict the analysis to the 
EPIC-pn data. The EPIC-pn camera was always operated 
in Small Window mode, with thin aluminum fllters to block 
visible light. The event files were reprocessed from the OPF 
data files using the epchain pipeline tasks of the XMM- 
Newton Science Analysis System (SAS version 6.5) and us- 
ing the most updated version of the public calibration files. 
These event files were then filtered for good time inter- 
vals following the "recipe" given in the XMM-A''ew;ton SAS 
handbook (V2.01 23 July 2004, Sect. 5.2.4). Taking into ac- 
count the dead time, the filtered event files then contained 
5.9, 7.6, 9.4, 30.0 and 12.9 ks of good exposure in the pn 
detector for OBS 1, 2, 3, 4 and 5 respectively. The final 
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Fig. 1. Top: the 0.5-10 keV X-ray light curve of the different EPIC-pn observations of Mkn 841. It varies by a factor ^ 
5 in 4 years. Middle: the 3-10 keV X-ray light curve of the different EPIC-pn observations of Mkn 841. It varies by a 
factor two times lower than the 0.5-10 keV light curve meaning that the broad band variability is dominated by the soft 
(< 3 keV) band one. On the other hand, flux variability up to ~50% is also observed on tens of ks in the soft and hard 
X-rays. Bottom: Hardness ratio light curve (5-10 keV)/(3-5 keV). 



EPIC-pn count rates for each observation are reported in 
Tab.[TJ They are always well below the 1% pile- up threshold 
for both instruments. Due to the proximity of the "Small 
Window" edges, the source spectra and light curves were 
built from photons detected within a 40 arcsec extraction 
window centered on the source. X-ray events correspond- 
ing to pattern < 4 were selected. The background was esti- 
mated within a window of the same size as the source from 
an offset position. 

In the following, all errors refer to 90% confidence level 
for 1 interesting parameter (Ax^=2.7). 



3. Model-independent analysis 

3.1. Light curves and hardness ratios 

We have plotted in the upper and middle panel of Fig.[T]the 
0.5-10 keV and 3-10 keV EPIC-pn count rate light curves 
of the different XMM-Newton observations of Mkn 841 as 
well as the hardness ratio (5-10 keV)/(3-5 keV) in the lower 
panel. The time binning is 500 sec. The total 0.5-10 keV 
count rate decreases by a factor ~5 in 4 years while the 
hardness ratio increases, reaching maximum values during 
OBS 4. Between 2001 and 2005, the 3-10 keV count rate 
shows variations a factor two lower than the total count 
rate meaning that the broad band count rate variability is 
dominated by the soft (< 3 keV) X-ray variability, at least 
on long time scale. Smooth soft and hard flux variabilities 
up to ^^50% are also visible on timescale of tens of ks. 



3.2. RMS 

Figure [2] shows the RMS spectra of the different observa- 
tions. The RMS func tion has been calcu lated following the 
procedure detailed in lPonti et al.l (|2004l ). We use the 0.5-10 
keV energy band. In order to increase the statistics we have 
grouped OBS 1, 2 and 3 to produce a single "2001" RMS 
spectrum. The time binning is of about 1 ks for each vari- 
ability spectrum while the energy binning has been chosen 
in order to have negligible Poisson noise. 
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Fig. 2. RMS spectra of OBS 1/2/3, OBS 4 and OBS 5. 



Due to the low degree of variability and the statistics we 
were not able to produce a spectrum with a large number 
of energy bins especially at high energy. Nevertheless the 
different RMS spectra are relatively flat: OBS 1/2/3 and 
OBS 5 are consistent with a constant value of ^^4 and 5% 
respectively. On the other hand the RMS spectrum of OBS 
4 is inconsistent with a constant (at more than 99.99%) 
with a mean value of 12%. In this case we note a rough 
increase of the RMS from 0.5 to about 2-3 keV and then a 
clear decreases at higher energies down to a few percent. 

From now on, due to the flux and spectral variability 
observed in OBS 4, we divide this observation in 3 parts 
(noted parti, part2 and part3) of about 15 ks duration each 
as indicated on Fig. [TJ Each part is analyzed separately. 

3.3. Spectral ratios 

We have plotted in Fig.[3]the ratios of the different EPIC-pn 
spectra with the EPIC-pn spectrum of OBSl. The binning 
is such that each bin has a 5 ct confidence level. No strong 
spectral variations are observed on short hours) time 
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Fig. 3. Ratios between the different EPIC-pn spectra and 
the EPlC-pn spectrum of OBSl. The dashed line corre- 
spond to the position of the 6.4 keV hne in the source frame. 

scale between OBS 1, OBS 2 and OBS 3 or the three parts 
of OBS 4. However, the 3-8 keV spectral shape clearly hard- 
ens between 2001 and 2005, in agreement with the increase 
of the hardness ratio shown in Fig. [TJ Above 8 keV the spec- 
tra keep roughly constant at least in shape. This spectral 
variability is more apparent in Fig. [4] where we have plotted 
the ratio between OBS 1 and the complete data set of OBS 
4. For comparison, we have over-plotted in this figure the 
expected variability assuming two different cases: a power 
law pivoting around 8 keV whose photon index decreases 
by 0.4 (dot-dashed line) and a variable black body peaking 
at 0.14 keV and varying by a factor 3 in flux (dashed line). 
These values have been chosen just to illustrate what kind 
of variability we can expect in such cases. Clearly neither 
of these simulations is able to reproduce the complete spec- 
tral variability suggesting a more complex behavior. We can 
also notice in Fig. [3] some variations of the iron line complex 
on long (i.e. year like between OBS 1 and part 3 of OBS 
4) and short (i.e. hours like between the different parts of 
OBS 4) time scale. 

4. Spectral analysis 

4.1. Phenomenological analysis 

Our first step in the spectral analysis was to fit the data 
with a simple power law, excluding the data below 3 keV. 
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Fig. 4. Ratio between the EPIC-pn spectrum of the com- 
plete OBS 4 and the EPIC-pn spectrum of OBSl. The data 
(white grey lines) have been rebinned in order to have a 5 tr 
confidence level or at least 25 data points per bin. We have 
over-plotted in this figure, the expected variability in two 
different cases: a pivoting power law around 8 keV whose 
photon index decreases by 0.4 (blue, dot-dashed line) and 
a variable black body peaking at 0.14 keV and varying by 
a factor 3 in fiux (red, dashed line). 
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Fig. 5. Ratio data/model for part 1 of OBS 4. The model 
is a simple power law fitted between 3 and 10 keV and 
extrapolated down to low energies. A strong soft excess 
and a line near 6 keV are clearly apparent. Note also the 
presence of a narrow feature near 4.8 keV. It will be discuss 
in Sect. O 



Figure [5] shows the ratio between the XMM-A^ew;ton/EPIC- 
pn data of part 1 of OBS 4 with the corresponding best fit 
power law. We clearly observed a strong soft excess below 
2 keV and a fluorescent iron line complex near 6.4 keV. 

The second step was to include very simple spectral 
components to reproduce the observed features. We use a 
power law for the continuum and a gaussian for the iron 
line. A more precise analysis of the line complex is done in 
the next section. We fit the data above 3 keV first. Then we 
fix the different parameters and include the data below 3 
keV down to 0.5 keV and we add a simple multicolor accre- 
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Table 2. Best fit values obtained with a simple model including a power law for the continuum, a gaussian for the iron line 
and a multicolor disc for the soft excess. The gaussian energy is given in the source frame. The hydrogen column density, 
not shown in this table, is always consistent with the galactic one, i.e. 2.34x 10^° cm~^. The Ax^ values correspond to 
the change in when adding the gaussian component. 



tion disc component (discbb in xspec) to model the soft 
excess. The best fit parameter values obtained with this 
method for the continuum, the multicolor disc and the line 
are reported in Table [S] 

The best fits are clearly never satisfactory and large dis- 
crepancies are present especially below 3 keV due to the bad 
black body approximation for the soft excess. Nevertheless 
several remarks, weakly affected by the goodness of the fit, 
can still been made: 



— The photon index reaches values as small as 1.3 during 
OBS 4 where the flux is the lowest. 

— The flux variation in the soft band (< 3 keV) between 
2001 and Jan 2005 is more than a factor 2 larger than 
the flux variation above 3 keV . 

— The iron line width appears variable on relatively short 
time scale (< 15 ks) e.g. between part 1 and part 2 of 
OBS 4. 



The spectral variability agrees with the results shown previ- 
ously with modcl-indepcndcnt methods and suggest a spec- 
tral pivot at high energy. However the presence of a sim- 
ple power law continuum is unlikely given the unusually 
(for a Seyfert galaxy) hard spectral index of OBS 4. More 
physical models are discussed in Sect. 14.31 Concerning the 
iron line, the 2005 observations confi rm the apparent line 
variability already obse rved in 2001 ((Petrucci et al.l . l2002t 
iLonginotti et al.l l2004D . The best fits reported in Table [1] 
indicate the presence of a narrow line (compared to the 
EPIC-pn resolution) except OBS 3 and part 2 of OBS 4 
where a broad component is preferred. The contour plots 
(at 68 and 90 % confidence level) of the line flux vs. line 
width for the different observations are plotted in Fig. [S] for 
2001 (left plot) and 2005 (right plot). The line variability 
is discussed in more detail in the next section. 

To decrease the number of degr ees of freedom, we choos e 
to follow the fitting procedure of ILonginotti et al] (|2004f ). 
We fit together observations with roughly the same under- 
lying continuum shape (cf. Tab. [2]) , e.g. the three pointings 
of 2001, on the one hand, and the different part OBS 4, on 
the other, keeping the power law continuum constant in 
shape but not in flux between the different spectra. OBS 5 
is analyzed as a single observation. 



4.2. The iron line complex 
4.2.1. Rapid variability 

The rapid (~hours time scale) variability of the line com- 
plex seems to be a common characteristic of this source. 
For a better visualization of this variability we have plot- 
ted in Fig. [7] the excess map of the complete set of the 
XMM-iVewton/EPIC-pn observa tions. This map has been 
obtained following the method of Iwasa wa et al.l 12004") and 
iTombesi et al- (2007^ apphed for NGC 3516 and NGC 3783 
respectively. We use resolutions of 2 ks in time and 100 eV 
in energy. Each temporal slice represents the residuals ob- 
tained when fltting each 2 ks spectrum with a power law 
between 3.5 and 8 keV but ignoring data between 4 and 7 
keV. This map is useful to reveal narrow features but sig- 
natures of broad ones can be mixed up with the underlying 
continuum. A narrow feature close to 6.4 keV is clearly 
present during most of the different pointings, but it seems 
to disappear from time to time on a very short time scale 
(- ks) hke in OBS 3 or in the middle of OBS 4 and OBS 5 

To check if the narrow line component variability is real 
or not we produce the light curves of its flux and equivalent 
width. They are plotted in Fig. |5| with a 4 ks binning time 
scale. The model includes a power law and a narrow (cr = 
eV) gaussian line with energy fixed at 6.4 keV. We fit the 
data between 3 and 10 keV. From time to time the narrow 
line is poorly detected with only upper limit on its flux and 
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Fig. 8. Equivalent width (top, in eV) and flux (bottom, 
in 10^^ ph cm"^ s~^) light curves of a narrow {a=0 eV) 
gaussian line whose energy is fixed to 6.4 keV. We use a 4 
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Fig. 6. Contour plots (68 and 90 %) of the line width vs. line flux obtained for the 3 observations of Jan 2001 (left) and 
the 3 parts of Jan. 2005 as well as Jul. 2005 (right). The model includes a simple power law + gaussian line. 
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Fig. 7. Excess map on time-energy plane of the different XMM observations. The energy scale is in the lab frame. 



EW. However both are consistent with a constant (at more 
than 98%) from 2001 to 2005. The fact that the EW is also 
constant is not surprising given the relative constancy of the 
underlying continuum near 6.4 keV (see Fig. [3]). A constant 
narrow line flux suggests the presence of remote reflection 
and it is discussed in the next section. If this interpretation 
is correct, the apparent variability shown by the excess map 
may reveal changes in either the continuum or any broad 
line emission underlying rather than changes of the narrow 
line itself. This is discussed in Sect. 14.2.31 



4.2.2. Remote reflection? 

The presence of a roughly constant narrow line component, 
as shown in Fig. [SI suggests the presence of neutral reflec- 
tion from remote material. Such reflection may come from 
the outer part of the accretion disc or even farther away 
from e.g. the dusty torus surrounding the nucleus as ex- 
pected in the unification framework of AGNs. 

In Tab.[3]we report the best fit parameters of the narrow 
((7=0 eV) gaussian line added to the 3-10 keV power law 
best fit of the three observations of 2001 keeping the power 
law flux free to vary between OBS 1, 2 and 3. We apply 
the same method for the three parts of OBS 4. OBS 5 
is analyzed normally as a single observation. The contour 
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6.20 6.30 6.40 6.50 6.60 
Line energy (keV) 

Fig. 9. Contour plots (68, 90 and 99%) line flux vs. line 
energy (source frame) of the narrow (a =0 eV) line compo- 
nent added to mimic emission from remote material. OBS 
1/2/3: red/solid contour, OBS 4 part 1/2/3: blue/dot-dot- 
dot-dashed contours, OBS 5: green/dot-dashed contours. 
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1/2/3 

4 part 1/2/3 
5 
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"■•^^-0.10 
0-^^-0.06 
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Table 3. Best fit parameters of the narrow (cr =0 eV) 
gaussian line added to mimic emission from remote mate- 
rial. The line energy is given in the source frame and the 
line flux in 10~^ ph cm~^ s^^. 



plotted the ratio of the complete OBS 4 data spectrum 
and the best flt power law -I- diskline -|- narro w (cr=0 eV) 
gauss ian line model. The diskline model (F abian et al.l 
Il989l) is a relativistic accretion disc line model around a 
static (Schwarzschild) black hole whose parameters are the 
line energy, the inner and outer disc radii, the power law 
index of the disc emissivity law and the disc inclination. 

This ratio has been obtained by fitting the data above 3 
keV and fixing the inclination to 30 degrees. This figure is 
then obtained after setting the line normalizations to zero. 
The presence of a complex line profile is clearly visible with 
a broad component down to ~ 4.5 keV. The solid line over- 
plotted in this figure corresponds to the diskline -t- narrow 
line shape. 

The corresponding best fit values are a power law pho- 
ton index T = 1.52l;gQ5, a diskline energy EpeK^ — 
6.12t°;^^ keV, a fine EW of 280ii^^g eV, an inner disc ra- 
dius Tin ~ 10.2'!i^i^ Tg and a disc emissivity power law 
index q > 2.4. The accretion disc outer radius is fixed to 
1000 Tg. The fit is very good with a x^/dof=147/165, sug- 
gesting that relativistic effects could indeed be a good ex- 
planation of the line profile. We apply this model to the 
different observations of Mkr 841. During the simultane- 
ous fit of OBS 1, OBS 2 and OBS 3 not only T but also q 
ad rin are kept constant between the different observations. 
We let the power law normalization as well as the diskline 
line flux and energy free to vary. We apply the same pro- 
cedure for the three parts of OBS 4. Concerning OBS 5, 
it is analyzed normally as a single observation. We fix the 
inclination angle to 30 deg and the outer radius to 1000 Vg. 
We obtain good fits in all cases. The corresponding best fit 
parameter values are reported in Tab. ID 



plots of the narrow line fiux vs. line energy are plotted in 
Fig. [HI The line fluxes obtained in the different observation 
periods are consistent with each other, in agreement with 
a roughly constant remote reflection. Surprisingly enough, 
while the best fit energy of the line is consistent with a 
slightly ionized iron line in 2005 (the 90% confidence level 
energy range being in between 6.43 and 6.47 keV for OBS 
4 and in between 6.43 and 6.54 keV for OBS 5) it is smaller 
than and marginally consistent with 6.4 keV (at less than 
10% confidence level) in 2001. 

In this last case, the detected line may be the signature 
of a slightly redshifted narrow iron line. This interpretation 
will be discussed in more detail in Sect. O If it is correct, 
the presence of remote reflection should still be tested for. 
We check the presence of a second narrow line component 
fixing the parameters of the first line to their best fit val- 
ues. The best fit parameters for this second narrow line are 
£;,,„e=6.48l°:J2keV and i^H„e=5.7t|^x 10-6 ph cm'^ s-\ 
These parameters are now in better agreement with those 
obtained for OBS 4 and OBS 5. 



4.2.3. Relativistic effects? 

If we believe in the constancy of the narrow line flux, the 
variability shown by the map excess (Fig. [7]) may result 
from variations in some other underlying component. The 
most common explanation of such variability is the presence 
of a broad and rapidly variable line component. This is also 
illustrated in the upper panel of Fig. (TU] where we have 
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Fig. 10. Ratio data/model for OBS 4. Top: The model is 
a simple power law -|- diskline -I- narrow ((7=0 eV) gaus- 
sian line fitted between 3 and 10 keV. We ignore the best 
fit line components (over-plotted in solid line) to produce 
this ratio. Bottom: The model is a simple power law -|- 
a partially ionized absorbing material with large velocity 
shear -I- narrow {a—0 eV) gaussian line. We also ignore the 
best fit line component (over-plotted in solid line) to pro- 
duce this ratio. The presence of a broad line component is 
strongly reduced and only the narrow component is present 
(see Sect. 112]). 
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Table 4. Best fits parameter values obtained witli a power law + relativistic line (diskline) + narrow {a—O eV) gaussian 
line model, fitting the data above 3 keV. OBS 1, 2 and 3 are fitted simultaneously keeping the power law index F, the 
disc emissivity power law index q and the inner disc radius ri„ constant between the different observations. We apply 
the same procedure for part 1, 2 and 3 of OBS 4. The line energy is given in the source frame, r^n is in unit of Vg and is 
larger than 6 (Schwarzschild metric). The inclination angle is fixed to 30 deg and the outer radius to 1000 Vg. 
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Table 5. Best fit parameter values obtained with a power law -f relativistic smeared absorption (swind) + narrow (ct=0 
eV) gaussian line model, fitting the data above 3 keV. OBS 1, 2 and 3 are fitted simultaneously keeping the power law 
index F constant between the different observations. We apply the same procedure for part 1, 2 and 3 of OBS 4. The 
absorbing column density Nh is in unit of lO^^cm^^. 



An alternative to the broad line component could be 
the presence of ionized warm absorption, smeared by rel- 
ativistic effects, whose imprints on the underlying contin- 
uum could mimic a broad line component (e.g. Reeves et al.l 
|2004|) . Such a model has also been suggested recently in the 
literature to produce the soft excess in AGNs and will be 
discussed more precisely in the next section, but we test 
this model here to see its impact on the iron line profile. 
We fit the 3-10 keV data with a simple power law and a par- 
tially ionized absorbing material with l arge velocity shear 
(swiN D model origina lly proposed by 'Gierlih ski fc Doni 
(120041 ) and updated bv lGierliiiski fc Do ne (2006;)). We also 
add a narrow gaussian line to mimic the remote reflec- 
tion. For comparison with the diskline model, we have 
plotted in the lower panel of Fig. [TO] the ratio data/model 
for the complete OBS 4 data spectrum. The presence of a 
broad line component is strongly reduced and only a nar- 
row component is observed. The corresponding best fit val- 
ues are a power law photon index F — I.OIqJ, a column 
density and ionization parameter of the absorbing material 
TV/j > 36 X lO^^cm-^ and ^ = 308012?to and a gaussian 
velocity smearing a — 0.19to o5- The fit is statistically ac- 
ceptable with a x^/dof=176/168. While it is larger than the 
DISKLINE fit, both are statistically equivalent following the 
F-test . We apply this model to the different XMM-Newton 
observations. We obtain good fits in all cases. The corre- 
sponding best fit parameter values are reported in Tab. [51 

In conclusion, the data are consistent with the presence 
of a constant narrow line, potential signature of a remote re- 
flection. Then the apparent line variability shown by the ex- 
cess map (Fig. [7]) may result more likely from slight changes 
of the underlying continuum. We have tested two possible 
origins of such variability. It could be due to the presence of 
a rapidly variable broad iron line component. It could also 
result from variable relativistically smeared absorption fea- 



tures. Both cases require strong relativistic effects to agree 
with the data. 



4.3. A broad band physical analysis 

Up to now, our spectral analysis was relatively phenomeno- 
logical, using simple components to fit the different spectral 
features present in the data (soft excess, broad and narrow 
line) and focusing on the high (> 3 keV) energy range. The 
next step is the use of more physical and consistent models 
on the total energy range of the EPIC-pn instrument. 

Recent studies suggest that relativistically blurred (pho- 
toionized) reflection from the accretion disc could be an 
appealing explanatio n for the presence of strong soft ex- 
cesses in AGNs (e.g. ICrummv et al.l |2Q06|) . Moreover, in 
the cases in which a broad Fe line is clearly detected 
(such as in MCG-6-30-15) the model is very robust be- 
cause the soft excess and broad Fe line are fitted self- 
consist ently with the same rel ativistically blurred reflection 
model. ICrummv et al.l (|2006f ) already applied this model 
to OBS 3 with success. As said in the previous section, 
an alternative explanation for the origin of the soft ex- 
cess assumes the presence of absorption features. Since 
the soft excess does not show strong edges nor absorption 
lines, strong velocity gradients are need ed in the absorbing 
medi u m to smear out the s e features (iGierliiiski fc Done! 



20061 : 



'2OOI; 'Sobolewska fc Don4 [20071: ISchurch fc Pond . 
Gierlihski & Done, 2006). A disc wind could produce such 
spectral signatures. However, this model does not directly 
explain the presence of a broad line feature and, if present, 
it has to be produced by another component. On the other 
hand, as shown in the previous section, ionized absorption 
may have some impacts on the continuum spectral shape 
near 6 keV and consequently on the observed line broadness 
(cf. Fig.[Tr 
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Our XMM-Newton data of Mkr 841 appear to be well 
adapted to test both interpretations. 

4.3.1. The fitting procedure 

Our base line model has the following components: 1) a 
neutral absorption free to vary above the galactic value, 2) 
a cut-off power law continuum (the high energy cut-off be- 
ing fixed to 300 keV) and 3) a neutral reflection, if needed, 
to reproduce the narrow line component. For the neutral 
reflection we use the tables of the Ross & Fabian code 
oss fc Fabianl . 120051 ) and we fix the ionization parameter 
to 1 to account for a neutral medium and the illuminating 
power law continuum to 1.9. This value corresponds to the 
aver age 2-10 ke V X-ray photon index of Seyfert galaxies 
(e.g. lMattl[200ll ). but the fit results do not depend signifi- 
cantly on this parameter. The free parameters of the first 
three components of our model are the hydrogen column 
density, the cut-off continuum power law photon index and 
normalization and the normalization of the neutral reflec- 
tion component. 

Then we add a fourth component to reproduce the soft 
excess i.e. either a relativistically blurred ionized reflection 
or a relativistically smeared ionized absorption (noted and 
ABS respectively in the following). For the ionized reflec- 
tion, we also u se the tables of the Ross & Fabian code 
oss fc Fabianl . 2005). The blurring is done by convolv- 
ing the reflection spectrum with a Laor profile (kdblur 
kernel). The Laor model is a relativistic accretion di sc line 
mode l around a maximally rotating (Kerr) black hole (jLaoii 
Il991h . The corresponding free parameters of are the inner 
radius r™ of the reflecting accretion disc, the disc emissiv- 
ity power law index q, the normalization and the ionization 
parameter ^ref of the ionized reflection component. The 
iron abundance is fixed to the solar one for the computation 
of the reflection, the disc outer radius to 1000 Schwarzschild 
radii and the inclination angle to 30 degrees. 

Concer ning ABS we use again t he SWIND model de- 
veloped by iGierhhski fc Pond ()2006[) for XSPEC. The free 
parameters are the column density N^^^ and ionization 
parameter ^abs of the absorbing wind and the gaussian 
velocity smearing (tabs- 

We apply these different models to the 0.5-10 keV en- 
ergy range of the EPIC-pn detectors. Like in the previous 
section, we analyzed OBS 1, 2, 3 simultaneously. We keep 
all the parameters constant between the different data sets 
except the normalization of the cut-off power law contin- 
uum and, for REF the normalization and ionization param- 
eters of the ionized reflection, and, for ABS, the column 
density and ionization parameter of the absorbing wind as 
well as the gaussian velocity smearing. We apply the same 
procedure for the three parts of OBS 4. Then we have 14 
free parameters for and 15 for ABS. OBS 5 is analyzed as 
a single observation with 8 and 7 free parameters for and 
ABS respectively. 

4.3.2. The best fit results 

While these models reproduce relatively well the global 
spectral shape of the data, the fits are not very good. The 
residuals show some narrow features especially in the soft 
band and suggest the presence of absorption. These features 
are even more pronounced in 2005. Such absorbing material 



is indeed confirmed by a quick look to the RGS data with 
the presence of absorption features due to the Unresolved 
Transition Array of Fe VII-XII. There is also apparently 
some variability between OBS 4 and OBS 5 indicating some 
variations of the absorber properties. A detailed analysis of 
this component is out of the scope of the present paper and 
will be done in a forthcoming publication (Longinotti et al. 
in preparation). For the present analysis, we model these 
absorption features by adding a warm absorber (WA) com- 
ponent (absori in xSPEC) in our fits letting the absorber 
hydrogen column and ionization state free to 

vary. The improvement is highly significant. For the Ax^ is 
30, 72 and 75 for the addition of two parameters for OBS 
1/2/3, OBS 4 part 1/2/3 and OBS 5 respectively For ABS 
the Ax^ is equal to 15, 106 and 70. 

The best fits for both models become now statistically 
acceptable. The corresponding best fit parameter values 
for the power law continuum, the neutral refiection and the 
WA are reported in Tabs [5] and [7] for and ABS respectively. 
The best fit parameter values characterizing the ionized re- 
flecting or absorbing material are reported in Tabs. [5] and [H]. 

Interestingly, not only the best fits with and ABS are 
statistically acceptable, they are also statistically similar 
and cannot be ruled out with the present XMM-Newton 
data. However they predict different spectral characteris- 
tics. For instance, the continuum spectral variability be- 
tween 2001 and 2005 is significantly larger with (Ar=0.9) 
compared to ABS (Ar=0.3). The former case is relatively 
unusual for a Seyfert 1 galaxy and would require strong 
changes in the emitting regions. The properties of the WA 
are also different, at least for OBS 1/2/3 and OBS 4, the 
column density and ionization parameters being larger by a 
factor 10 (and even more for ^'^^) between and ABS. On 
the other hand, the neutral reflection normalization is con- 
sistent with a constant between the different pointings with 
both models, supporting the presence of remote reflection. 

For comparison we have plotted the unfolded best flts 
obtained with and ABS on Fig. [TT] for OBS 1, the flrst 
part of OBS 4 and OBS 5. We have also over- plotted the 
different spectral components like the neutral reffection 
and, for the model, the relativistically blurred ionized 
reflection. Large differences above 10 keV are expected 
especially for 2001 and Jan. 2005. 

We can also analyze more precisely the results ob- 
tained for each model separately, beginning with REF. 
Interestingly, while this model requires still extreme values 
of the disc inner radius, the disc emissivity power law index, 
now also constrained by the soft excess, is close to 4 in OBS 
4 part 1/2/3 and OBS 5, in better agreement w ith theoret- 
ical expectations (e.g. 'Mart occhia et all |2000[ ) . This con- 
trasts with the results obtained with the DISKLINE model 
where larger values of q are generally found (see Tab. [4|). 
We agree that these discrepancies can be partly explained 
by intrinsic differences between the Laor and diskline 
profiles. Note however that q is still very large (> 8.6) in 
OBS 1/2/3 whatever the relativistic line profile used. The 
photon index is also very steep (F ^^2.45) in these observa- 
tions. Noticeably, constraining q to be smaller than 5 gives 
a more reasonable photon index with F =2.261qq5 while 
the other parameters do not significantly change. The is 
larger (xVdof=804/787 i.e. Ax^=20) but still acceptable. 
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Table 6. Best fit parameter values for the continuum, the neutral reflection and the WA obtained with the relativistically 
blurred ionized reflection model. is the neutral hydrogen column density, F the continuum photon index, Nneutrai the 
neutral reflection normalization, the ionized hydrogen column density of the WA and its ionization parameter 

.Ax^ is the fit improvement with the addition of the warm absorber. The corresponding best fit parameter values of the 
reflecting material are reported in Tab. [51 
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Table 7. Best fit parameter values for the continuum, the neutral reflection and the WA obtained with the relativistically 
smeared ionized absorption ABS model. The parameters definition is the same as in Tab.[6l The corresponding best fit 
parameter values of the absorbing material are reported in Tab. [HI 
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Table 8. Best fit parameter values characterizing the ionized reflecting material in the rhodel. ^ref is the ionization 
parameter of the ionized reflection component, q the disc emissivity power law index and the inner radius of the 
reflecting accretion disc. The reflection fraction is the ratio of the ionized reflection flux in the 0.1-1000 keV band divided 
by the total 0.1-1000 keV flux of the best flt model. 



We have also reported in Tab. [8] the reflection fraction, 
i.e. the ratio between the flux in the ionized reflection 
component and the total flux. It is of the order of 30-50% 
which is relatively large since we expect ^ 10-20 % in the 
case of an isotropic illumination. Finally the properties of 
the WA significantly change for OBS 5 compared to the 
other observations. 



On the other hand, the characteristics of the WA are 
relatively constant between the different pointings with 
the ABS model. The continuum power law index is also 
close to the standard one for a Seyfert galaxy i.e. '^1.9, 
reaching also a steeper value 2.2) in 2001. To flt the 
data the ABS model requires a change in the proper- 
ties of the relativistically smeared absorbing material 
between 2001 and 2005, the column density and the veloc- 
ity smearing being larger in 2001. These parameters are 
more consistent with each other between OBS 4 and OBS 5. 



4.3.3. Comparison with a simultaneous BeppoSAX 
observation 

The BeppoSAX instruments cover the 0.1-200 keV range 
and thus better constrains the reflection component com- 



pared to XMM-Newton. Mkr 841 was pointed by BeppoSAX 
between the 11th and 14th of January 2001, i.e. partly si- 
multaneously with the XMM-Newton observations of 2001, 
for a total net time exposure of ~90 ks for the MEGS 
instrument (2-10 keV energy range) ~40 ks for the PDS 
(10-200 keV energy range) and only 20 ks for the LEGS 
(0.1-2 keV energy range). Since the BeppoSAX observation 
w as much longer than the XMM-Newton one (see Fig. 1 
of iPetrucci et al.ll2002f) . we did not try to fit the XMM- 
Newton and BeppoSAX data altogether. We only checked 
the consistency of the BeppoSAX data above 10 keV with 
the XMM-Newton best fit model expectation. 

The unfolded best fit spectra obtained with and ABS 
for OBS 1 are plotted in Fig. [H] with the BeppoSAX 
MEGS and PDS data (we ignore the LEGS data due to 
their poor quality). Only the MEGS normalization was 
let free to vary, the PDS one being fi xed to 0.86 the 
MEGS one (following iFiore et all |1999[ ). The agreement 
is relatively good given the possible variability of the 
reflection component during the BeppoSAX pointing and 
the fixed value of 300 keV of the power law continuum high 
energy cut-off in the fits of the XMM-Newton data. We 
note however a better agreement of the ABS model with 
the BeppoSAX /PDS data above 10 keV. With the rhodel, 
the flux above 10 keV is underestimated by a factor ~2. 
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Table 9. Best fit parameter values characterizing the absorbing material in the ABS model. N^^^ and S,abs are the 
column density and ionization parameter of the absorbing wind and aABS the gaussian velocity smearing. 
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Fig. 11. Comparison of the unfolded best fits ob- 
tained with the and ABS models for OBS 1, the first 
part of OBS 4 and OBS 5. The thick sohd fines cor- 
respond to the total best fit model and the dashed 
thick line to the ABS one. Over-plotted in these fig- 
ures in thin lines are the different components of each 
model. Thin sohd lines: cut-off power law continuum, 
relativisticaUy blurred ionized reflection and neutral 
reflection. Thin dashed lines: cut-off power law con- 
tinuum, modified by relativisticaUy smeared ionized 
absorption, and neutral reflection. The fluxes are in 
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However this depends on some of our model assumptions. 
For example constraining q to be smaller than 5 instead 
of being completely free (see Sect. I4.3.2p gives a better 
agreement (the green data points on the left plot) with 
the BeppoSAX /PDS data, similar to the agreement found 
with ABS. 



5. Indication of redshifted narrow iron lines 

5.1. A strongly redshifted line at 4.8 keV in OBS 4 

We focus here on a possible detection of a narrow line at 
4.8 keV in the first part of OBS 4. This feature is clearly 
visible on Fig. [S] but also on the excess map shown on Fig. 
[71 To estimate its confidence level we first added a gaussian 
with energy 4.8 keV to a simple power law best fit model 
of OBS 4 part 1. The gaussian energy and flux are let free 
to vary but, given the apparent narrowness of the feature, 
we flx the gaussian width to zero. The excess is clearly 



detected with an improvement of the fit of Ax^ = 11. The 
best fit parameters of the line are Eune = 4.80+q q3 keV 
(source frame) and an equivalent width EW = 50 ± 20 eV. 
The contour plot of the line energy versus the line flux is 
plotted in solid line on Fig. [T31 

We also estimate the confldence level of this detection 
by simulating a large number of spectra with a continuum 
similar to that of OBS 4 part 1. We flt the different 
simulations with a simple power law. We then add a 
gaussian, flxing its energy between 4 and 8 keV by step of 
0.1 keV and flt again the data, the gaussian normalization 
being free to vary while the gaussian width was flxed 
to 0. Then for each simulated spectra, we keep the best 
Ax^ resulting from this procedure. For 1000 simulations 
we only flnd 14 cases where a Ax^ >11 corresponding to 
a confldence level of 98.5%. The strong detection in the 
EPIC-pn data is however attenuated by the weak detection 
in the MOS data. Indeed we have plotted in dashed line 
on Fig. [S] the corresponding contour plot when fltting the 
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Fig. 12. Comparison of the OBS 1 best fit obtained with the (left) and ABS (right) models with the simultaneous 
BeppoSAX. data. The different spectral components are over- plotted in the upper plots: left: neutral reflection (dashed 
line), ionized blurred reflection (dotted line), cut-off power law (dot-dashed line) and total (solid line). The grey (green 
on the colored version) ratios correspond to the best fit ihodel with the disc emissivity power law index q forced to be 
smaller than 5 (see Sect. 14.3^ : right: neutral reflection (dashed line), cut-off power law modified by smeared absorption 
(dot-dashed line) and total (solid line). 



EPIC-pn and MOS data simultaneously. The detection is 
now significant at only 84.4% from our simulations. 

If we assume that this line is a redshifted neutral fluo- 
rescent iron line, its observation at 4.8 keV (source frame) 
implies a redshift factor of ~ 0.75. Such high redshift value 
suggests an origin close to the central engine where rel- 
ativistic effects become important. Moreover, the narrow 
shape of the line constrains the emitting region to be suffi- 
ciently small. There are growing evidences of such variable 
narrow emission featur es in the X-r ay spectr a of severa l 
AGN in the literature (iTurner et all. [20 04: Por guet et"al] . 
12004 iPetrucci et al.l . [20021: iTombesi etal.1 . 120071. The com- 
mon scheme to explain these components suppose the pres- 
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Fig. 13. Contour plot energy (source frame) vs flux of the 
line detected close to 4.8 keV during part 1 of OBS 4. The 
contours correspond to 68, 90 and 99% The solid and dash- 
dotted contours correspond to the EPIC-pn and EPIC-pn 
-|- EPIC-mos data respectively 



ence of transient magnetic flares briefly illuminating a lo- 
calized part (hot spot) of the accretion disc and producing 
the iron line by fluorescence. 

Following Longinotti et al. (2004, Fig. 4) we have pro- 
duced a 90% probability map (from the 90% error on the 
line energy) of the accretion disc where the line could come 
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Fig. 14. A schematic map of the inner accretion disc show- 
ing the allowed regions where a narrow line with parame- 
ters consistent with the one observed during part 1 {Eune — 
4-79to o2 keV) could be emitted. In red: a =0 and i =25. 
In blue: a =0.998 and i =25. In cyan: a = 0.998 and i =60 



deg. The black circle correspond to a radius of 6 r^. The 
disc is rotating counterclockwise. 
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Fig. 15. Contour plot (68 and 90 %) energy (source frame) 
vs flux of the narrow line ((7=0 eV) added to the DISKLINE 
fit of OBS 1 (sohd red hne), OBS 2 (dotdotdot-dashed blue 
line) and OBS 3 (dot-dashed green line). The narrow line 
energy in OBS 1 peaked close to 6.2 keV (source frame) 



from, for different values of the black hole specific momen- 
tum and disc inclination angles (cf. Fig. I14p . For an incli- 
nation of 25 deg. the line emission can come from region at 
a distance as small as 3 rg from the central black hole and 
even closer for larger inclinations. Moreover, if we assume 
that we observed the same hot spot on the disc for more 
than an orbit then the global line profile is expected to be 
relatively complex (e.g. iDovciak et aLll2004[ ) with a strong 
blue peak. If we interpret the observed feature wi th this 
blue peak, a Kerr solution is unavoidable (|Pechacek et al.l . 
I2005D . 

5.2. A slightly redshifted line at 6.2 keV in OBS 1 

We discuss now the possible detection of a narrow line fea- 
ture near 6.2 keV (source frame) in the first observation 
of 2001 as already noted in Sect. 14.2.21 Starting from the 
best fit DISKLINE model reported in Tab. 21 we add a nar- 
row ((T=0 eV) gaussian line letting its energy free to vary 
for each data set OBS 1, OB S 2 an d OBS 3. This differs 
from the analysis done in Sect. l4.2^ where the narrow gaus- 
sian line parameters were kept constant between the differ- 
ent observations. The corresponding energy vs. flux contour 
plot of each narrow line are reported in Fig. [121 The nar- 
row line energy is consistent with the neutral fluorescent 
iron line energy of 6.4 keV in OBS 2 and 3 but is incon- 
sistent at more than 90% with this value for OBS 1 where 
Eune = 6.23t°:°^ keV. 

If interpreted, as above, as a (slightly) redshifted nar- 
row iron line we can also produced a 90% probability maps 
of the accretion disc where the line could come from, for 
different values of the black hole specific momentum and 
disc inclination angles. This has been done in Fig. 1161 

6. Summary and Discussion 

We have presented in this paper a detailed spectral analysis 
of the XMM-Newton/EPlG-pn data of Mkr 841 including 
all the XMM-Newton observations of this object. Given the 



20 r 




Fig. 16. Same as Fig. [T31 but for a narrow line with pa- 
rameters consistent with the one observed during OBS 1 
{Eiine = 6.23j^Q Q3 keV). We have limited ourselves to the 
inner part (< 20 rg) of the accretion disc. Cf. Fig.[T31for the 
color scheme. The blue and red regions overlapped almost 
totally. 



relative spectral complexity of this source we would like to 
summarize first the different results we obtained. 

Mkr 841 has been observed three times in January 2001 
(OBS 1, OBS 2 and OBS 3) and two times in 2005 in 
January (OBS 4) and July (OBS 5). Flux and spectral 
variations being present during the ~ 45 ks of OBS 4, this 
observation has been divided in 3 parts for the spectral 
analysis. A strong soft excess as well as a complex iron line 
proflle, both known to be present in this source for a long 
time, are clearly detected in all these pointings, the high 
sensitivity of the XMM-A''ewfon/EPIC-pn unveiling very 
puzzling spectral and temporal behavior. The 0.5-10 keV 
flux varies by a factor 3 in 4 years and is dominated by 
the soft band (< 3 keV) variability, the data above ~5 keV 
keeping roughly constant on short and long time scales. 
The spectral variability is also important. Fitting the 3-10 
keV data with a simple power law the spectral index varies 
from ~1.9 in 2001 down to ^1.3 in Jan. 2005. However a 
simple pivoting power law component cannot explain by it- 
self the broadband 0.5-10 keV spectral variability observed 
on year time scale thus indicating a more complex spectral 
variability. 

The line profile is also complex, being apparently a 
mixture of broad and narrow components. It is highly 
variable and the 2005 data confirrn the rapi d line vari- 
ability observed in 20 01 by iPetrucci et all (|2002D and 
iLonginotti et al.l ()2004l ) with a variability time scale as 
short as a few kilo-seconds. This strong X-ray variability 
suggests a small emitting region close to the black hole. 
Fitting the line with a diskline model gives a good repre- 
sentation of the line profile and requires in some cases steep 
disc emissivity power law indexes. 
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Given the high complexity of the spectra, we choose 
to decrease the degree of freedom of some of our fits by 
fitting simultaneously observations with roughly the same 
underlying continuum, e.g. the three pointings of 2001 as 
well as the different parts of OBS 4. In these cases, during 
the fits the power law continuum was kept constant in shape 
but not in flux between the data sets. Concerning OBS 5, 
it was analyzed normally as a single observation. A narrow 
line component is present in all the pointings and its flux 
is consistent with a constant between 2001 and 2005. This 
suggests the presence of remote reflection in the data. If 
this interpretation is correct, then the observed variability 
of the iron line complex is more likely due to the variability 
of the underlying component below the narrow iron line, 
either due to a broad line component or strong smeared 
absorption. 

Then we analyze the broad band (0.5-10 keV) EPIC- 
pn energy range. We applied two different models, recently 
proposed in the literature, for the origin of the soft ex- 
cess: a relativistically blurred ionized reflection (REF) and 
a relativistically smeared ionized absorption (ABS)). We 
also added a neutral reflection to reproduce the narrow line 
component. Both models are a reasonable representation of 
the overall broad band shape but are formally statistically 
unacceptable du to the presence of absorption features in 
the soft energy range. The addition of a warm absorber 
strongly improves the fits that converge to statistically ac- 
ceptable and statistically equivalent representations of the 
data. Both models are also consistent with the partly si- 
multaneous BeppoSAX observations done in 2001. 

Finally we also note the presence of a strong narrow 
feature near 4.8 keV in the first part of OBS 4. Its 
detection is marginally significant in the EPIC-pn data 
(>98.5%) but the significance decreases to 84% when we 
include the MOS data. If interpreted as the blue horn of a 
relativistically distorted neutral iron line, the large redshift 
implies the presence of a Kerr black hole. 

The most remarkable result of our study is to reproduce 
reasonably well the broad band (0.5-10 keV) and complex 
spectral characteristics of Mkr 841 with a small number 
of spectral components. However we were not able to dis- 
criminate between the two different interpretations of the 
soft X-ray excess, which have been debated in the recent 
literature. Moreover both models are able to reproduce the 
spectral shape close to the iron line either by adding a broad 
component, like in the case of the model, or by modifying 
the power law shape, like in the the ABS model. This is ex- 
emplified in Fig. [IT] where we show a zoom of the unfolded 
best fit models obtained with and ABS for OBS 5 in the 
1-8 keV range. 

Whatever the model used, the need of strong relativis- 
tic effects to smear/blur the absorption/emission features 
suggest an origin of the X-ray emission close to the central 
black hole. Then the absorbing material required in the 
ABS model could be the base of a disc wind/jet surround- 
ing the inner X-ray emitting region. The external part of 
this wind, with lower velocities, may then explain the ab- 
sorption features of the WA. 

On the other hand, the model is based on the presence 
of strong refiection as indicated by the large refiection 
fraction that we obtained to fit the data. Another way of 
quantifying this large amount of reflection is to measure 
the commonly used reflection parameter R. Since this is 
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Fig. 17. Unfolded best fits of OBS 5 obtained with (solid 
line) and ABS (dashed line) in the 1-8 keV range. The 
curvature of the spectral shape below to the iron line, due 
to a broad line component in and by smeared absorption in 
ABS, is present in both cases. 



not a direct parameter of the refiection tables that we used, 
we estimate its value by roughly adjusting the reflection 
shape (ignoring the blurring effects) with the one expected 
with the PEXRiv model of xspec fixing the power law 
photon index and normalization as well as the ionization 
parameter of pexriv to our best fit values. Following this 
method, we obtain R values of the order of 10, 2 and 4 for 
OBS 1/2/3, OBS 4 and OBS 5 respectively. We recaU here 
that the blurring effects, as well as the ionization level of 
the refiecting material considerably smooth the reflection 
shape. In consequence even a reflection component with 
i? = 10 is strongly attenuated in the outgoing spectrum. 
But such large reflections necessarily require the presence 
of some effects that increase the reflection component 

compared to the illuminating co ntinuum like e.g. in 

the light bending effect s model (jMiniutti et al.l 120031 : 
iMiniutti fc Fabianl . l2004j) or the i nhomogeneous accre- 
tion flows model proposed recently bv lMerloni et al.l (|2006[ ). 

The low statistics of our RMS spectra (cf. Fig. ^ do not 
allow a detailed comparison with the variability behavior 
expected with ABS or and thus unable us to discriminate 
the two models. We can note however that, while both mod- 
els c an correctly reproduce the RMS spectra of some AGNs 
fe.g. lGierlinski &: Donell2006l : iPonti et al.ll2006H . ABS pre- 
dicts an RMS spectrum which rapidly increases and peaks 
aro und 1 keV, where t he ab sorber enhances the variabil- 
ity (iGierlihski fc Pond . [2006[ ) . This is not seen in our RMS 
spectra and so argues against this model a bit or at least 
that there is no absorber variability on short timescales. 
But this has to be tested on data of better quality. 

Finally, we note that our (admittedly marginal) detec- 
tion of redshifted narrow iron lines in the data better fits 
in the relativistically blurred reflection interpretation where 
such features are naturally expected. But we agree that this 
is a relatively weak argument that requires further investi- 
gations. 
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7. Conclusion 

Mkr 841 is a bright Seyfert 1 galaxy which is known to 
possess a strong soft excess and complex iron line profile. 
We have presented a detailed analysis of the whole XMM- 
Newton/EPlC-pn data of this object from 2001 to 2005. 
The conclusions of this analysis can be summed up as fol- 
lows: 

— Strong flux and spectral variability are observed on 
month and year time scales. The 0.5-10 keV flux varies 
by a factor 3 in 4 years and is dominated by the soft 
band (< 3 keV) variability, the data above ^ 5 keV 
keeping roughly constant on short and long time scales. 

— The iron line is apparently a mixture of broad and nar- 
row components. Its profile is rapidly varying on very 
short time scale (a few ks). 

~ The broad band 0.5-10 keV spectrum is well described 
by a model including 1) a neutral absorption, 2) a cut- 
off power law continuum, 3) a neutral reflection and 
a fourth component for the soft excess. We use two 
different models for this last component: a relativisti- 
cally blurred photoionized reflection (ihodel) and a rel- 
ativistically smeared ionized absorption (ABS model). 
Both models give statistically acceptable and statisti- 
cally equivalent fits of the data even including, for 2001 
observations, partly simultaneous BeppoSAX. data up to 
200 keV. 

— Both models agree with the presence of remote reflec- 
tion characterized by a constant narrow component in 
the data. However they differ on the presence of a broad 
line component present in but not needed in ABS. 
Consequently, the physical interpretation of the line pro- 
file variability is quite different, resulting from the vari- 
ability of the broad line component in and from the 
variability of the absorbing medium in ABS). 

— Mkr 841 is also a good candidate for the observation 
of redshifted narrow iron lines, two marginal detections 
being discussed in this paper. If such features would 
favor the relativistically blurred photoionized reflection 
interpretation, their marginal detections do not permit 
any clear conclusion. 

While it seems reasonable that the reality is a com- 
plex combination of absorption and reflection effects the 
present analysis underlines the difficulty in well disentan- 
gling these components in the X-ray spectra of AGN. Broad 
band observations including data above 10 keV (Suzaku, 
INTEGRAL) will certainly help to better test both inter- 
pretations but the high sensitivity in the 10-50 keV en- 
ergy range is crucial. We may need to wait for missions like 
SIMBOL-X or XEUS to have a clear understanding of what 
processes really take place in sources like Mkr 841. 
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